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ABSTRACT 

HII regions surrounding supermassive black holes (BHs) in an otherwise still neutral 
intergalactic medium (IGM) are likely to be the most easily detectable sources by 
future 21cm experiments like LOFAR. We have made predictions for the size distribu- 
tion of such HII regions for several physically motivated models for BH growth at high 
redshift and compared this to the expected LOFAR sensitivity to these sources. The 
number of potentially detectable HII regions does not only depend on the ionisation 
state of the intergalactic medium and the decoupling of the spin temperature of the 
neutral hydrogen from the cosmic microwave background (CMB) temperature, but is 
also strongly sensitive to the rate of growth of BHs at high redshift. If the supermas- 
sive BHs at redshift 6 were built up via continuous Eddington-limited accretion from 
low mass seed BHs at high redshift, then LOFAR is not expected to detect isolated 
QSO HII regions at redshifts much larger than 6, and only if the IGM is still signifi- 
cantly neutral. If the high-redshift growth of BHs starts with massive seed BHs and 
is driven by short-lived accretion events following the merging of BH hosting galax- 
ies then the detection of HII regions surrounding supermassive BHs may extend to 
redshifts as large as 8 — 9 but is still very sensitive to the redshift to which the IGM 
remains significantly neutral. The most optimistic predictions are for a model where 
the supermassive BHs at z > 6 have grown slowly. HII regions around supermassive 
BHs may then be detected to significantly larger redshifts. 
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1 INTRODUCTION 

Reionisation - the change of state of the intergalactic 
medium (IGM) from a neutral state at recombination to the 
present day ionised state - is an epoch which is currently un- 
der intense study. Analyses of the spectra of high redshift 
QSOs indicate a smooth increase in the optical depth of the 
IGM with redshift which accelerates at 2: > 6.1 (Fan et al. 
2002, Songaila 2004, Fan et al. 2006). Quasars at 2 > 6.1 
show almost complete absorption of the flux blue-ward of 
the Lyman— a resonance [the Gunn-Peterson (1965) effect]. 
Inferring the ionisation state of the IGM at these high red- 
shifts from QSO absorption spectra is difficult and ambigu- 
ous (e.g. Oh & Furlanetto 2005, Bolton & Haehnelt 2006). 
If the small amounts of flux transmitted just blue-ward of 
the Lyman— Q centre are attributed to the presence of a 
locally ionised region embedded within a significantly neu- 
tral medium, then lower limits on the neutral fraction of 
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iHii may be inferred (Wyithe & Loeb 2004a; Mesinger & 
Haiman 2004). Malhotra & Rhoads (2005), however, argue 
that the number density of Lyman-a emitters can be used 
to place a lower limit on the ionised fraction of gas at 2 ~ 6 
at 20 — 50 per cent. Fan et al. (2006) also come to the con- 
clusion that the IGM probed by the highest redshift QSOs 
is predominantly ionised. 

The optical depth to Thompson scattering inferred from 
the recent WMAP measurements offers an independent con- 
straint on the integrated reionisation history. Interpreted 
as a sudden reionisation at a redshift 2:roion the year three 
WMAP results give 2rcion = lOtj.s (Spergel et al. 2006). 

The emerging picture is consistent with the reionisa- 
tion of cosmic hydrogen being complete by 2 ~ 6 — 13 and 
therefore significant contribution to reionisation by sources 
beyond current detection limits. Experiments to detect the 
redshifted 21cm signature of high redshift neutral hydrogen 
will potentially provide copious amounts of information on 
the epoch of reionisation. The European design for a low 
frequency interferometer, LOFAR, may be in operation as 
early as 2007 and other low-frequency experiments (see 
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are in various stages of planning and development. However, 
mapping the high redshift neutral gas distribution in 21cm 
will be extremely challenging. 

The redshifted 21cm signal due to ionised regions pro- 
duced by bright, high redshift QSOs surrounded by an oth- 
erwise still neutral intergalactic medium may be our best 
bet for an detection of individual objects due to their 21cm 
signature (see Tozzi et al. 2000 for a detailed discussion). 
Kohler et al. (2005) used numerical simulations to show that 
ionised regions around high redshift QSOs are indeed likely 
to be the most prominent features in the high redshift neu- 
tral gas distribution. Wyithe, Barnes & Loeb (2005) use a 
semi-analytic model for the high redshift QSO HIT region 
distribution and predict detectable HII regions at 2 ~ 6 — 8. 
However, the growth history of high redshift supermassive 
BHs is very uncertain. In this paper we will investigate a 
range of models for high redshift BH growth which repro- 
duce the bright end of the luminosity function of redshift 6 
QSOs and sample the range of plausible growth histories at 
higher redshift. 

We begin by reviewing the physics of detecting a 21cm 
signal due to neutral hydrogen (^. In section|3we describe 
our models for the evolution of the QSO luminosity func- 
tion (QLF). We discuss the expected size of HII regions due 
to QSOs in [21 Finally, in ^ we compare our computed 
size distributions to the expected sensitivity of LOFAR^ to 
extended sources and discuss observational prospects. 

Throughout this paper we adopt a cosmological matter 
density Qrn = 0.27, baryonic matter density fib = 0.044, cos- 
mological constant Qa = 0.73, present day Hubble constant 
Ho = lOO/i km s"^ Mpc"^ = 71km s"^ Mpc"\ a mass vari- 
ance on scales of 8h~^ Mpc as = 0.84 and a scale invariant 
primordial power spectrum (slope n — 1). 



2 PHYSICS OF 21CM 

EMISSION/ ABSORPTION 

A ground state neutral hydrogen atom can emit or absorb 
at a wavelength of 21cm corresponding to the energy of the 
transition between the singlet and triplet hyperfine levels. At 
the low frequencies of this transition, the specific intensity of 
the emission is proportional to the brightness temperature 
of the source Tj, as described by the the Raleigh- Jeans law, 

/(/.) = ^kBTb. (1) 

The brightness temperature of a typical patch of sky 
at redshift z depends on the optical depth of the gas to 
21cm radiation, the temperature of the cosmic microwave 
background Tcmb and the spin temperature of the gas, 

r6 = TcMBe-" + (l-e-")T,, (2) 

where r depends on the density and ionisation state of the 
gas. This results in a differential brightness temperature. 
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Significantly neutral hydrogen gas for which Ts differs 
from Tcmb will emit or absorb at 21cm in contrast to the 
CMB. The value of Ts is determined by the kinetic tempera- 
ture of the IGM, TiGM, and the strength with which Ts and 
TiGM are coupled via coUisional ionisation and/or Lyman— a 
pumping. TiGM is expected to decouple from Tcmb as early 
as 2 ~ 150 — 200 as the timescale for Compton scattering off 
residual free electrons becomes longer than the expansion 
timescale (e.g. Scott & Rees 1990). In the absence of a heat 
source, the IGM cools adiabatically below Tcmb and there- 
fore the very first (Population-Ill) sources are expected to 
be born in a cold IGM. 

The redshift range over which HII regions will appear 
in absorption is uncertain; if Ts couples to Tigm prior to sig- 
nificant long range X-ray heating there may be a significant 
absorption epoch. Furlanetto (2006) argues that if the (ther- 
mal) radiation from Pop-II stars dominate reionisation then 
this will be the case, but if Pop-Ill stars were dominant in 
reionising the universe then the IGM could be heated above 
Tcmb before there are sufficient Lyman— a photons to cou- 
ple Ts to Tigm (see also Chen & Miralda-Escude 2003 and 
Kuhlen, Madau & Montgomery 2006). For simplicity we as- 
sume here that the temperature of the IGM lies well above 
the temperature of the CMB. Regions of ionised gas then 
appear as holes in the 21cm emission due to the IGM (see 
also Kohler et al. (2005) for a more detailed discussion). 

With the rise of the galaxy and QSO populations the 
ionised fraction of gas increases and the 21cm signal emission 
due to the IGM becomes weaker. The signal will disappear 
by 2 < 6 when we know the neutral hydrogen is predom- 
inantly ionised. There will be a finite period of time when 
the IGM is heated above the temperature of the CMB but 
most of the gas is still neutral. This is the (uncertain) red- 
shift range over which instruments like LOFAR can probe. 
In addition, the performance of LOFAR will depend on our 
ability to model and remove foreground contamination due 
to the ionosphere, our own galaxy, and both extended and 
point-like extragalactic sources. 

3 MODELS FOR THE LUMINOSITY 

FUNCTION OF HIGH REDSHIFT QSOS 

3.1 The growth of supermassive black holes at 
high redshift 

We are interested here in the HII regions around supermas- 
sive BHs at high redshift. Accretion onto a BH is an efficient 
mechanism for producing ionising photons and therefore a 
supermassive BH that has been built up by gas accretion 
should generate a very large ionised bubble. If the ionis- 
ing flux generated by the supermassive BHs at high redshift 
exceeds that produced by other sources in the same region 
then the sites of supermassive BHs will show up prominently 
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in 21cm maps as large ionised regions in an otherwise still 
neutral IGM. In order to estimate the space density and size 
distribution of these HIl regions we require a model for the 
growth of the supermassive black holes and the production 
of ionising photons. Our best anchor point for such mod- 
elling is the observed luminosity function of z = 6 QSOs. 
We therefore consider here models which reproduce this lu- 
minosity function but differ in the timescale over which the 
black holes powering the z = & QSOs grow. 



3.2 Passive evolution model 

The e-folding time for BH growth at the Eddington limit 
is ~ 0.44(eacc/(l — Eacc)) Gyr, where eacc is the fraction 
of the mass energy released as radiation. Estimates for the 
mass of seed BHs produced by Population-Ill remnants are 
~ 10 — 10"^ M0. Such seed masses would require ~ 15 — 20 
e- foldings to grow to 3 x lO" Mq , the estimated mass of the 
black holes powering the observed z = G QSOs. The growth 
of a 3 X 10^ M© BH by a redshift ~ 6 from a PopuIation-III 
remnant therefore requires more or less continuous accretion 
if the accretion rate is limited to the Eddington accretion 
rate. It is a matter of intense debate whether accretion above 
the Eddington rate contributes significantly to the growth of 
supermassive black holes (see Begelman, Volonteri & Rees 
2006 for a recent discussion). 

We first construct a simple passive evolution model for 
the growth of supermassive black holes and the evolution 
of the QSO luminosity function at 2; > 6 in which a fixed 
comoving density of black holes accrete at their Eddington 
hmit with a fixed duty cycle /duty 

The bolometric luminosity of a QSO powered by a BH 
of mass Mbh radiating at the Eddington limit is given by 



hfetime of a QSO in this model is redshift dependent, and 
depends on the duty fraction of the luminous population, 



tq = flumtniz) 



3.5 X 10 yrs I 
for z » 1. 



l + z 



-3/2 



(7) 



3H(z) 



is the Hubble time at a redshift z and 



Here tniz) 

H{z) = Ho[{l + z)^nm + Oa]^''^ is the Hubble parameter. 
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3.3 A merger driven model 

3.3.1 Galaxy mergers as triggers for QSO activity 

The activity powering optically bright QSOs is thought to 
be triggered by galaxy mergers, which drive cold gas onto a 
central black hole. Various models for the growth of black 
holes in hierarchically merging galaxies have been built (e.g. 
Kauffmann & Haehnelt 2000, Volonteri et al. 2003) and have 
been successful in reproducing the observed QSO luminos- 
ity function at 2; < 5 . These models generally assign a QSO 
light curve to each merging event, guided by observation- 
ally and/or physically motivated estimates for the central 
black hole mass and associated accretion rate (see Hopkins 
et al. 2005 for light curves motivated by numerical simu- 
lations of a galaxy merger). Below we construct a simple 
version of such a model for the QLF which uses the (ex- 
tended) Press-Schechter formalism for dark matter halos 
to estimate the rate of galaxy mergers (see Erickcek et al. 
(2006) for a recent critical assessment of the accuracy of the 
extended Press-Schechter formalism). The model is similar 
to those described by Haehnelt, Natarajan & Rees (1998) 
and Wyithe & Loeb (2003). 
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where is the Thompson scattering cross section. The B- 
band luminosity is assumed to be a fraction /s = 1/10 of 
Lboi. Assuming that a fraction eacc = 0.1 of the mass energy 
of the accreted mass is released as radiation, the ensemble 
average luminosity evolution with respect to the redshift 6 
population is then given by. 
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5 X 10 yrs- 
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where tz = t{z) is the cosmic time at redshift z. 

We adopt the QLF of H3.3I as our template QLF at 
2; = 6. In this model all luminosities grow exponentially 
and the shape of the luminosity function is independent 
of redshift. We consider two values of the total duty cy- 
cle /duty ~ [0.1, 1.0]. To account for a possible contribution 
to the growth in an optically obscured accretion mode we 
split the duty cycle of accreting black holes into phases of 
luminous growth (/lum) and growth in an obscured phase 
(fobs). We assume the QSO spends equal amounts of time 
in each phase so that /lum ~ fobs ~ [0.05, 0.5]. The observed 



3.3.2 Estimating merger rates of black holes 

Our model makes a number of simplifying assumptions. We 
assume that the merger rate of black holes is equal to the 
merger rate of the dark matter halos of the galaxies hosting 
the black holes as estimated using the (extended) Press- 
Schechter formalism. We assume a relation of the mass of the 
BH, Mbh, and the mass of the dark matter halo forming in a 
merger M as motivated by the observed correlation between 
Mbh and the velocity dispersion of the host galaxies stellar 
bulge /dark matter halo. 



Mbh OC cr". 



(8) 



in nearby galaxies (Gebhardt et al. 2000, Ferrarese & Mer- 
ritt 2000, Ferrarese 2002). This relation appears to hold in- 
dependently of redshift for QSOs out to 2 ~ 3 (Shields et 
al. 2003, 2006). Empirical estimates of a typically fall in the 
range 4 — 5. We have chosen a = 5 which is consistent with 
a simple self-regulated growth scenario in which the black 
hole grows until it radiates enough energy to unbind the gas 
that is feeding it (e.g. Silk & Rees 1998, Haehnelt Natarajan 
& Rees 1998, Wyithe & Loeb 2003). We further assume that 
the velocity dispersion of the stellar bulge is proportional to 
the circular velocity of the dark matter halo at the virial 
radius, a = v^ir/V^, with (Barkana & Loeb 2001), 
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Figure 1. The thin curves show the QLF as a function of bolometric QSO luminosity at z = 6, 8, 10, 12, 20 for the merger rate model 
with Co = 10~^ * {left panel), and the passive evolution models with obscured duty fractions equal to the duty fraction of the luminous 
population and total duty fraction /duty = 0.1 (middle panel) and /duty = 1-0 {right panel). The thick line shows the Fan et al. (2004) 
fit to the observed z ~ 6 luminosity function at a rest frame wavelength of 1450Jl, where we have assumed that 10 per cent of the energy 
is radiated in the B-band and converted to a B-band magnitude using Mb = Afi45o — 0.48 + 5logio ^ ■ Note that the dashed line 
in the bottom left corner of panel 3 is the 2 = 8 QLF. 
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where Q,^ - n„(i+z)3+nA+nfc(i+z)2 ' 
IStt^ + 82d — 39d^ is the overdensity of a virialised halo at 
redshift z. Equations lO and lead to a dependence of 
black hole mass on host halo mass M and redshift of the 
form, 
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Note that to sets the ratio of Mhh/M for a given redshift. 

Assuming that all dark matter halos host a central BH 
with mass given by equation IIIO^ , the merging rate of halos 
of mass {M — AM) and AM is equal to the merger rate 
of halos hosting BHs of mass (Mbh — AMbh) with those 
hosting BHs of mass AMbh for AMbh = e{AM,z{t))AM 
and Mbh = e{M,z{t))M. A''haio(M, AM, t) is given by the 
product of the probability per unit time that a halo mass 
AM will merge with another halo to form a halo with mass 
M and the space density of halos with the appropriate mass 
difference. 



iVhalo(M, AM,t) 



dn 



dAMdt d{M - AM) 

iVbh (Mbh, AMbh, t). 



(11) 



The mass density is derived from the Press-Schechter (1974) 
halo mass function with the modification described by Sheth 
& Tormen (1999) and the halo merger probability is taken 
from the extended Press-Schechter formalism (Lacey & 
Cole 1993). 

3.3.3 From merger rate to luminosity function 

We assume that following a merger, each black hole accretes 
and radiates at its Eddington limit for a period equal to the 
dynamical time of gas in a disk of radius 0.035 times the 
virial radius as in Wyithe & Loeb 2003, 



= tdyn{z) 

= lO'' yrs 



-1/2 



1+Z 



'3/2 



(12) 



The light curve is assumed to be a simple top-hat function 
L(Mbh,t)/L0 = rMbh/MQe[tq-(t-tz)] (with 9 the usual 
Heaviside step function). 

The QLF can then be estimated by integrating over the 
merger rate of all BH pairs as follows. 



ML,z) = 



/■O.SeAf 

dMbh / dAMbh 
xN^^{M^^,AM^i„t')5[L - L(Mbh,t)]. 



^ ,dt' 

az 

dz' 



(13) 



If we perform the integral over Afbh using the 5 func- 
tion, use equation 1101 to make the change of variable from 
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BH masses to halo masses and assume that the merger rate 
changes Httle over the QSO hfetime [i.e. tq « tniz)], we 
get (see Wyithe & Loeb 2003), 

rO.SM o 

= / dAMN{M,AM,t{z))tq— (14) 

Jo.25M 

1 

£ L0 M-i ' 

where L/Ijq = Ct{M, z)M/}Aq. Note that the integral over 
AM is truncated at the lower end to ensure that we are only 
counting major mergers and at the upper end to count each 
potential halo merger pair only once. 

In principle, a, eo and tq may be tuned to fit our model 
QLF to the observed QLF at z ~ 6. We have fixed a = 5 
and tq — tdyn{z) and varied to to fit to the observed lumi- 
nosity function. For to = 10~^'^ we obtain an acceptable 
fit to the observed luminosity function (see Figure 1). This 
normalisation is consistent with the corresponding empiri- 
cal normalisation of equation (|SJ for galaxies in the local 
universe (eo ~ 10"" ■^ see Ferrarese 2002). 

Note that with the assumption tq = t^yn the QSO life- 
time at 2 = 6 is « 1.3 X 10^ yrs, which is consistent with 
other estimates of the duration of the optically bright phase 
of QSOs (Haehnelt, Natarajan & Rees 1998, Haiman & Hui 
2001, Martini & Weinberg 2001, Hopkins et al. 2005). Note 
that the shorter lifetime in the merger driven model com- 
pared to the passive evolution model results in smaller HIT 
regions for the assumed values of /duty. 



3.4 Comparison with observations 

3.4-1 Quasar luminosity function 

Figure 1 shows the evolution of our model QSO luminos- 
ity function at 2: > 6 (the passive evolution model with 
two values of the duty cycle and the merger driven model). 
By construction all models agree with the observed space 
density of bright QSOs at z ~ 6 as determined by Fan et 
al. (2004). The redshift evolution differs dramatically for 
the three models. Fan et al. (2001) find that the observed 
space density of bright QSOs decreases with redshift like 
g-i.i5z j-pj. 2 g < 2 < 6. We compare this to the redshift 
evolution in our model QLFs by summing over the bright 
end of each QLF. We sum from a minimum luminosity 
iboi = 5.3 X 10^'^ L0 which corresponds to a rest frame 
1450A absolute magnitude of —26.8 [for a conversion to B- 
band magnitude Mb — M1450 — OA8 + 5logioj^ (Fan et al. 
2001)]. The results are displayed in Figure 2. For the merger 
driven model the density of QSOs declines with a similar 
slope to the decline in the observed space density at lower 
redshift. In the passive evolution model the rate of decline 
depends sensitively on the assumed duty fraction. For the 
lower value of the assumed duty cycle (/duty = 0.1) the black 
holes grow relatively little at high redshift and the drop in 
the space density of bright QSOs is significantly slower than 
in the merger driven model. For the high value of the duty 
cycle (/duty ~ 1) the QSOs are growing rapidly and quickly 
decline in luminosity at higher redshift. As we will discuss 
in more detail below the three models span a wide range of 
possible growth histories of supermassive black holes at high 




z 

Figure 2. Redshift evolution of QSO space density for each of 
the three models plotted in Figure 1; the merger rate model with 
lifetime of tq = tdyn(z) from equation 1121 {thin solid line), the 
passive evolution model with fixed duty cycle /duty =0.1 {dashed 
line) and /duty = 1-0 {dot-dashed line). The solid thick line shows 
the evolution found for SDSS QSOs in the redshift range 3.6 < 
z < 5 (see Fan et al. 2004) and the dot is the estimated space 
density of QSOs at 2 = 6 (Fan et al. 2004). We have taken Lmin ~ 
5 X 10^^ Lq which is the bolometric luminosity corresponding 
M1450 = -26.8 assuming that Mb = M1450 -0.48 + 5Zogio(Q^) 
and ten per cent of the bolometric luminosity is radiated in the 
B-band. 

redshift. We should also note that all three models overpre- 
dict the space density of bright QSOs at lower redshifts. As 
discussed in detail by Bromley, Somerville & Fabian (2004) 
the rapid growth history of supermassive black holes at high 
redshift is not easily reconciled with the much slower ob- 
served growth at low redshift which is normally modelled 
with rather severe feedback effects due to energy input by 
stars and/or the black hole. As we are only interested in the 
growth of supermassive black holes before reionisation we 
make no attempt to model this here. 

3.4-2 Black hole mass densities and accretion rates 

The accretion rate of mass onto black holes during the op- 
tically bright phase of accretion can be extracted from the 
optical QLF as follows, 

f^Pacc,lum J / , / dt f dn , f \ T / /\ 1 tacc 

(15) 

An approximation to the total mass density in BHs can 
be extracted from the model luminosity function by dividing 
the luminosity function by the fraction of luminous/active 
halos at a given time /ium(i, z), 

C Lq Mq J dL /ium(i,2) 



6 K.J. Rhook & M.G. Haehnelt 



I 



N 
X) 



O 



fduty = 0-1 

1.0 




20 



Figure 3. Redshift evolution of tlie mass density in black holes 
{bold lines) and observed accretion rate of matter onto black holes 
{thin lines) as described in i|ij.4l for the three models described in 
El The hatched regions shows the range of estimates for the mass 
density in BHs at 2 = (Yu & Tremaine 2002). 



4 HII REGIONS AROUND HIGH REDSHIFT 
SUPERMASSIVE BLACK HOLES 

Assuming that the QSO is embedded in a homogeneous gas 
distribution with neutral hydrogen number density uh ~ 
3::hi'^h,o(1 + z)'^, the increase of the physical volume of 
ionised hydrogen Veil generated by an ionising flux Nionit) 
is given by (e.g. Madau, Haardt & Rees 1999, hereafter 
MHR99), 



dt 



3H{t)VHU = 



(17) 



We assume the cosmic abundance of hydrogen with a neu- 
tral fraction xm and that the temperature of the ionised 
IGM (Tigm) is 10'* K. The average hydrogen recombination 
timescale tree for such a gas is (see Abel et al. 1997), 



3 X 10 yrs 



/ TlGM 

I 10* K 



7 



Ccff 

3 



(18) 



where we have assumed an effective clumping factor CcB = 
3. The effective clumping factor relevant for recombinations 
is expected to be low since the densest self-shielded regions 
of neutral gas have a small covering factor and should not 
impede the progression of the ionisation front (see Miralda- 
Escude, Haehnelt & Rees 2000 for a detailed discussion). 

Equation 1171 is integrated numerically to solve for the 
physical radius of the HII regions for QSOs in each of our 
QLF models. 



For the passive evolution model /lum is assumed to be 
constant. Note, however, that this will somewhat overesti- 
mate the BH mass density. For the merger rate model fium 
is estimated to be the ratio of the QSO lifetime to the aver- 
age survival time of a dark matter halo as a function of halo 
mass and redshift [ts{M, z)] calculated from the probability 
distribution in equation (2.21) of Lacey & Cole (1993). At 
low redshift ts is approximately the Hubble time, but it de- 
creases to 30 per cent of the Hubble time by z ~ 12 for the 
relevant halo masses. Figure 3 shows the evolution of the 
black hole mass density and the rate of mass accretion dur- 
ing optically bright phases for our three models. As expected 
the evolution differs strongly between the three models. We 
also show the observed black hole mass density at z — 
which is about 30 — 300 times larger than our estimates at 
z = 6. The difference in the integrated observed accretion 
rate and the mass density in BHs is due to the assumed 
periods of obscured QSO accretion. By construction the in- 
tegrated luminous accretion rate will fall short of the mass 
density in black holes in the passive evolution models by a 
factor /duty //lum = 2 at 2 = 6. In our merger driven model 
this difference is much larger. Approximately ten per cent 
of the inferred BH mass density is accumulated due to lumi- 
nous accretion. In this model approximately 90 per cent of 
the growth would have to occur in an obscured phase. Such 
a growth scenario is consistent with the findings of Hopkins 
et al. (2005) based on hydrodynamical simulations of QSO 
activity fuelled by merging disk galaxies. 



Equation (1171 is approximately equivalent to assuming that 
the volume of the ionised region is equal to the ratio of 
the number of emitted ionising photons and the density of 
neutral hydrogen, 

Vii ^ ^ ■ (20) 

riH 

For a QSO with constant luminosity, this approximation 
breaks down as tq —> tree and the HII region reaches its 
Stromgren radius. However for a QSO with luminosity grow- 
ing exponentially with the mass of the central BH, most of 
the HII region growth occurs between t{z) — K and t{z). Since 
K ~ tree, the result holds for the QSOs with tq >t rcc in our 
models. Note that since little ionising flux is produced early 
in the lifetime of a QSO powered by a BH continually ac- 
creting at the Eddington limit, the size of the HII regions 
at 2 = 6 are approximately equal for each of the passive 
evolution models even though the lifetimes differ by a factor 
of 10. 

For each of our models the radial extent of the ionised 
region generated by a QSO depends inversely on the sur- 
rounding neutral hydrogen fraction xm and increases with 
the QSO luminosity as -Rhii oc L^/^x^j^'^ . 

We calculate Mon at all redshifts by assuming the op- 
tical/UV spectral energy distribution [I/(i^)] as given by 
MHR99 with a high energy cut-off of 200 keV. The spec- 
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tral energy distribution is normalised by assuming ubLb = 
fsLhoi (with fs = 0.1 as in This gives an ionising flux, 



iVion = 



^du ^ 1.5 X 10" s-^ 
hv 



10"L, 



(21) 



where i/ion is the minimum photoionising frequency. The cor- 
responding radius of the HII region, for a QSO with constant 
luminosity is, 



3.6 Mpc(2:Hi)"^/^ 

ta 



Lhoi 



1/3 



(22) 



-| 1/3 / , I s -1 
' 1 + Z 



1 X lO^rs 



5 PROSPECTS FOR DETECTING HIGH 
REDSHIFT HII REGIONS 

5.1 Future high-redshift 21cm experiments 

The next generation of radio telescopes are aiming to ex- 
plore the low-frequency universe, with detection of an epoch 
of reionisation signal a major scientific goal. Various low- 
frequency interferometers are being designed, tested and 
built with some already taking data. The largest radio tele- 
scope on the horizon is the Square Kilometre Array (SKA^). 
The SKA is still in the design phase and not expected to 
be in operation for more than a decade. However there 
are several smaller-scale telescopes including LOFAR (to 
be built in the Netherlands), the Mileura Wide-Field Ar- 
ray (MWA sited in Western Australia) and the Primeval 
Structure Formation telescope (PAST ^, currently being 
built in China) that are expected to be in operation on much 
shorter time-scales. These arrays are typically composed of 
~ few X 10'^ dipoles (with approximately 10 per cent of 
the collecting area and therefore sensitivity of the tentative 
SKA design) spread over a region with maximum baseline 
~ 1 km and designed to operate at subsets of the frequency 
range 50 — 300 MHz. Below we focus on the possibility of de- 
tecting HH regions with LOFAR, but our results are easily 
transferred to alternate designs for a low-frequency interfer- 
ometer. 



5.2 LOFAR sensitivity to HII regions 

We assume that the ionised region has a brightness temper- 
ature contrast of (JThii = 27 mK [(l-|-z)/10]^'^^, correspond- 
ing to the limit where the hydrogen in the surrounding IGM 
is completely neutral, and the spin temperature (Ts) of the 
surrounding neutral gas is efficiently coupled to the kinetic 
temperature of a warm IGM. 

The compact LOFAR core, rather than the full array, 
will offer the highest sensitivity to extended structures. In 
the following we have adopted array parameters correspond- 
ing to the published design study of LOFAR ^ . Assuming a 



^ http://skateIescope.org/ 

^ http:/ /haystack. mit.edu/arrays/MWA 

^ http://web.phys.cmu.edu/past/ 

^ http://www.lofar.org 




Figure 4. Signal to noise ratio as a function of physical radius 
for spherical HII regions at 2 = 6, 8, 10, 12, 20. The calculation 
assumes constant sky noise at 250 K, 100 hours integration time, 
a spectral resolution of 0.1 MHz, a top hat beam with the same 
diameter as the source, and an observed bandwidth which encom- 
passes the source. 



spectral resolution Avch and integration time Ar, an ar- 
ray of A'^stn dual polarisation antennae, each with collecting 
area AcS and system temperature Tsys, will have a bright- 
ness temperature sensitivity 



AT(,(_Bmax) = 



2fcBr 



128 mK 



250 K 



0.1 MHz 



-1/2 



Ar 
100 hr 



^max 


2 


( iVstn \ 


-1 




2 km 




V3200 ) 




18 m2 



where Bmax is the maximum distance between antennae in 
the core of the array, F = 57bcam/ (A/Bmax)^ and f^bcam is the 
solid angle of the synthesised beam, fcs denotes the Boltz- 
marm constant. The sensitivity to point sources o-p is given 
by, 



8.8 X 10"^ Jy 



250 K / VO.l MHz 



-1/2 



/ A. ^ 


-1/2 






V 100 hr ) 




18 m2 


V3200 ) 



LOFAR is expected to have a bandwidth of 4 MHz with a 
spectral resolution of 4 kHz which will likely be binned into 
~ 0.1 MHz channels. The spatial resolution of the LOFAR 
core will be inversely proportional to the maximum distance 
between the antennae (Sres ~ A/-Bmax ~ 2.5 arcmins^iy^). 
In comparison, the angular size (Omi) and frequency depth 
(Ai^Hii) of an HII region with radius -Rhii are given by 
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Rhii C^P^) ^hii C^P^) ^hii (^^P'^) 

Figure 5. The black lines show the size distribution of HII regions (physical radius) around active QSOs for z = 6,8,10,12,20 for 
the QLFs in Figure 1. The grey lines show the size distribution if fossil HII regions are included. The horizontal line corresponds to a 
source density of one per 400 deg'^ per 4 MHz of observed bandwidth. The vertical line corresponds to the smallest HII region at z = 6 
detectable in 100 hours of integration time, as calculated in >I5.2I however we find that the minimum detectable physical size changes 
very little with redshift. Note that the surrounding IGM has been assumed to be neutral. 



fen 



2^(1 + ,), 



14 arcmins 



Al/HII 



Rmi \ 
5 Mpc J 

H{z) Rmi 



l + z 



Dc(z) 
Dc{z = 6) 



(24) 



4.6 MHz 



f Run \ 




\ H(z) 1 


\ 5 Mpc ) 




H{z = 6)_ 



,(25) 



where Vo ~ 1420 MHz is the rest-frequency of the signal and 
Dc{z) is the comoving distance to a source at redshift z. 

We now calculate the signal to noise ratio (SNR) with 
which a spherical bubble of a given radius at a certain red- 
shift can be detected with A'^stn ~ 3200 antennae each with 
Aefi = 18 in a circle of variable diameter -Bmax. The 
brightness temperature sensitivity increases with the com- 
pactness of the array. The surface brightness sensitivity will 
be maximal when the size of the beam is the same as the 
angular size of the source -Bmax ~ A2icm/fen, F = 1. 

The SNR for a spherical HH region is then given by. 



PLC 



5Th 



1/2 



TYai. (26) 



max — A21cm/fen) V A/^ch / 

The factor Tyui accounts for the fraction of the beam that 
is filled with ionised gas. For a spherical top heat beam 

.^Hii ~ ^{^'L^"" = 2/3. Tyoi may be adjusted for different 
array patterns/synthesised beams and different HH region 
geometries. 

This prescription allows us to estimate the signal to 
noise ratio for detection of an HH region as a function of 
its redshift and radius. Figure 4 shows the calculated signal 
to noise ratio as a function of the physical radius of an Hll 



region at various redshifts. Note that the increase of SNR 
with size of the HH region will become shallower once the 
size exceeds the beam size for the most compact configu- 
ration possible. This calculation indicates that the smallest 
bubble surrounded by a neutral IGM that may be detected 
at the 5 — (T level is ~ 2.3 Mpc at « = 6. The minimum ra- 
dius changes very little with source redshift, and turns over 
at 2 ~ 10 due to the dependence on the angular diameter 
distance. 

Unfortunately, instrumental noise is not the only type 
of noise relevant for the detection of HII regions. Residu- 
als from the subtraction of extragalactic and galactic fore- 
grounds and the correction of the distortion of the signal by 
the ionosphere will all contribute to the noise. We currently 
have little idea how large the noise due to these residuals 
will be. The intensity fluctuations in the 21cm emission due 
to density fluctuations of the IGM will also act as noise (e.g. 
Kohler et al. 2005, Mellema et al. 2006). For 100 hour inte- 
grations with the LOFAR core the fluctuations due to the 
IGM surrounding the source limit the sensitivity on scales 
larger than ~ 6 arcminutes (Wyithe, Loeb & Barnes 2005); 
increasing the exposure time will not significantly decrease 
the minimum size of detectable HII regions. Furthermore, 
the edge of an HII region surrounding a QSO is not ex- 
pected to be sharp (see Zaroubi & Silk 2005) and the intrin- 
sic shape will appear distorted due to hght travel time effects 
(see Wyithe & Loeb 2004b; Yu 2005). The signal to noise 
ratios estimated here should thus be considered as (probably 
rather optimistic) upper limits. 



5.3 The space density of HII regions 

In Figure 5 we show our predictions for the space density 
of HII regions as a function of size at z > 6 in units of 
deg~^MHz~^ for each of the three models. The vertical and 
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horizontal lines show the typical size of a LOFAR observ- 
ing volume (~ 400 deg^ and 4 MHz of bandwidth) and our 
estimate for the 5 — a detection threshold for the size. The 
effective recombination timescale for hydrogen gas at 10* K 
is comparable to or somewhat shorter than the Hubble time. 
As discussed by Wyithe, Loeb & Barnes (2005) HII regions 
are thus expected to remain as 'fossils' for some time after 
the QSO turns off. We assume that the HH regions con- 
tinue to expand with the Hubble flow once the quasar turns 
off and there are no further recombinations, and count the 
number of fossil regions with physical radius R at z by sum- 
ming over the size distribution function at redshift intervals 
corresponding to the quasar lifetime. We assume that the 
quasars at 2 = 6 are on average half way through their op- 
tically luminous lifetime. There are essentially no fossil HH 
regions in the case of /duty = 1.0 since the QSO luminosity 
grows exponentially on a timescale comparable to the quasar 
lifetime. Note that for the other two models l//fossii(^, z) is 
~ 4 times smaller than l//ium for the smallest HH regions 
(corresponding to the faintest quasars, Lboi ^ 10^^ Lq) at 
2 = 6 due to the rapid evolution of the space density of 
QSOs. l//fossii(i?, z) asymptotes to one for large HH regions 
and high redshifts as expected. We show both the number of 
HII regions around active QSOs (black lines) and the total 
number (including fossil bubbles) of HII regions (grey lines) . 
We have assumed that the IGM is still neutral all the way 
down to 2: = 6. We will come back to this point below. 

For the merger-driven model we predict ~ 14 detectable 
sources around active QSOs per LOFAR field at 2: = 6. This 
number increases to ~ 340 if we include fossil HII regions 
around supermassive black holes which are no longer ac- 
tive. The numbers are similar to those reported by Wyithe, 
Loeb & Barnes (2005) for their similar model. In our merger 
driven model at least one detectable HII region per LOFAR 
field is predicted out to z ~ 7. Due to the steep slope of the 
QLF at high luminosities, the numbers would decrease more 
rapidly with redshift if the effective minimum detectable ra- 
dius is larger than we have estimated. 

For our two passive evolution models the predictions are 
quite different. The model with the lower duty cycle offers a 
more optimistic picture than the merger driven model. The 
longer QSO lifetime results in larger HII regions for the same 
QSO luminosity. We predict ~ 135 detectable HII regions 
per LOFAR field at 2 ~ 6 around active luminous QSOs and 
detectable HII regions are present up to 2 ~ 10. The value of 
the maximum redshift depends sensitively on the duty cycle 
and could be - at least in principle - much larger for smaller 
values of the duty cycle. This would, however, require the 
formation of very massive black holes at very high redshift 
which would be difficult in most formation scenarios. 

The passive evolution model with the higher value of 
the duty cycle depicts a much more pessimistic scenario. 
The predicted numbers at 2 = 6 are similar to the other 
models, but the the presence of detectable HII regions does 
not extend much beyond 2 = 6. This is easily explained by 
the very recent growth of the 2 = 6 supermassive black holes 
in this model. 

Figure 6 shows the redshift evolution of detectable ac- 
tive and fossil HII regions with 7?min = 2.3 Mpc. Clearly 
the prospects for detecting HII regions depend sensitively 
on how BHs grow at high redshift. 

So far we have assumed that the IGM is still neutral 



all the way down to 2 = 6. However the neutral fraction 
at redshift 6 is uncertain and the Universe may already be 
highly ionised at 2 ~ 6 — 7 (Malhotra & Rhoads 2005, Fan 
et al. 2006). A smaller but still significant neutral fraction 
would result in a larger HII region [see equation (I22ll 1 , how- 
ever the signal to noise ratio for the detection of a resolved 
HII regions scales with the square root of the neutral frac- 
tion (plc oc xl^j^); sources in a partially neutral IGM are 
more difficult to detect. 

For consistency we have calculated the volume filling 
factor of hydrogen ionised by the QSOs in each of our mod- 
els. We assume that our model extends to QSOs as faint as 
Lhoi = 10^^ Lq, include the contribution due to fossil re- 
gions and neglect recombinations. Our filling factor should 
therefore represent an upper limit on the expected value. 
We find a filling factor of HII regions generated by the UV 
emission of active quasars in our models which is less than 
~ 35 per cent at 2 ~ 6 for the model with the low duty 
fraction and for the merger model and ~ 5 cent for the 
model with the large duty fraction. For each model the fill- 
ing factor drops rapidly at higher redshift. Note again that 
we have integrated down to rather faint quasars and that 
taking recombinations into account would also lower these 
values. Our assumption that the HII regions generated by 
QSOs are isolated may therefore be somewhat questionable 
at 2 ~ 6 but should be reasonable at 2 > 7. Not surpris- 
ingly the quasars in our models which are consistent with 
the rather low observed space density of QSOs at 2 ~ 6 
alone are barely able to fully reionise hydrogen by 2 ~ 6. 

The quasars should also contribute to the ionisation 
of hydrogen with their X-ray emission. Due to their long 
mean free path. X-rays will lead to a more uniform ionisa- 
tion which includes the low-density IGM far away from the 
quasars (e.g. Madau et al. 2004, Ricotti & Ostriker 2004). 
The efficiency with which an X-ray photon ionises hydrogen 
atoms in the IGM depends on the average number of ioni- 
sations induced by secondary electrons, which is a function 
of the initial photon energy as well as the existing ionised 
gas fraction (see ShuU & van Steenberg 1985). We follow Ri- 
cotti & Ostriker and assume that each photon with energy 
> 1 keV results in 10 ionisations in a partially ionised IGM. 
Assuming again that our model extends to quasars as faint 
as Lboi = 10^^ Lq the mean ionised fraction at 2 = 6 due to 
the X-ray emission of active QSOs in our model is less than 
1 percent for all models. 

Another possibly relevant effect is the clustering of 
QSOs. Bright sources are sufficiently rare that clustering 
should not be important. For fainter sources overlap of clus- 
tered HII region could increase the sizes compared to our 
estimates in particular if fossil HII regions are taken into 
account. 



6 CONCLUSIONS 

We have investigated the prospects of using the 21cm sig- 
nature of the regions around supermassive black holes at 
2 > 6 to probe the growth history of these supermassive 
black holes and the ionisation history of their surrounding 
IGM. 

We have presented three simple, physically motivated 
models for the high redshift growth history of supermassive 
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Figure 6. Number of HII regions larger than the ealculated op- 
timal detection limit of iJmin = 2.3 Mpc detectable limit as 
a function of redshift for the merger driven model {solid line) 
and passive evolution models with /duty = 0.1 [dashed line) and 
/duty = 1 {dot-dashed line). The black lines correspond to HII re- 
gions around active QSOs and the grey lines include fossil HII re- 
gions. The horizontal line corresponds to one source per 400 deg^ 
field of view per 4 MHz bandpass. Note that the surrounding IGM 
has been assumed to be neutral. 

black holes in which the growth is either simple Edding- 
ton limited accretion or due to short-lived accretion phases 
triggered by the merger of galaxies. We have used a sim- 
ple analytical model for the size of HII regions in order to 
generate predictions for the size distribution of HII regions 
for each of our model QLFs. Our results show that mod- 
els which are all consistent with the luminosity function of 
bright QSOs at 2 = 6 but differ in the growth history of the 
supermassive black holes give very different results for the 
predicted number of detectable HII regions. 

Our estimates for the number of detectable HII regions 
at z ~ 6 are promising, up to several hundred detectable 
sources per LOFAR field around active QSOs and possi- 
bly many more fossil regions. However, the number of de- 
tectable HII regions depends sensitively on unknown details 
of SMBH growth and the ionisation state of the IGM. If 
the QLF evolves to z > 6 as predicted in our merger-driven 
model, then isolated HII regions will be prominent features 
of the high redshift neutral hydrogen distribution provid- 
ing reionisation was completed late {z < 7). However if the 
SMBH population is growing rapidly due to continuous ac- 
cretion, or reionisation occurred much earlier than 2; ~ 6, 
then the detection of QSO HII regions may prove elusive. 

The most optimistic picture emerges if supermassive 
black holes grow more slowly over an extended period of time 
as in our Eddington limited accretion model with a smaller 
duty cycle (/duty ~ 0.1). In this case detectable sources may 
be present to very large redshifts and the prospects for de- 
tection are good even if reionisation is completed relatively 
early. 



Observation of future 21cm experiments like LOFAR 
should not only provide a wealth of information on how the 
Universe was reionised but may also provide important clues 
for the growth history of supermassive black holes at very 
high redshift. 
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